Relativistic Hartree-Bogoliubov description 

of the deformed ground-state proton emitters 

D. Vretenar 1 ' 2 , G.A. Lalazissis 2 , and P. Ring 2 

1 Physics Department, Faculty of Science, University of Zagreb, Croatia 
2 Physik- Department der Technischen Universitdt Miinchen, D-85748 Garching, Germany 

00 , 

Q\ ' Ground-state properties of deformed proton-rich odd-Z nuclei in the re- 

gion 59 < Z < 69 are described in the framework of Relativistic Hartree 
Bogoliubov (RHB) theory. One-proton separation energies and ground-state 
quadrupole deformations that result from fully self-consistent microscopic cal- 
culations are compared with available experimental data. The model predicts 

>. 

£^ the location of the proton drip-line, the properties of proton emitters beyond 

m ' 

the drip-line, and provides information about the deformed single-particle 
qq \ orbitals occupied by the odd valence proton. 
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The decay by direct proton emission provides the opportunity to study the structure of 
systems beyond the drip-line. The phenomenon of ground-state proton radioactivity is de- 
termined by a delicate interplay between the Coulomb and centrifugal terms of the effective 
potential. While low-Z nuclei lying beyond the proton drip-line exist only as short lived res- 
onances, the relatively high potential energy barrier enables the observation of ground-state 
proton emission from medium-heavy and heavy nuclei. At the drip-lines proton emission 
competes with (3 + decay; for heavy nuclei also fission or a decay can be favored. The pro- 
ton drip-line has been fully mapped up to Z = 21, and possibly for odd-Z nuclei up to 
In [ffll. No examples of ground-state proton emission have been discovered below Z = 50. 
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Proton radioactivity has been studied from odd-Z nuclei in the two spherical regions from 
51 < Z < 55 and 69 < Z < 83. The systematics of proton decay spectroscopic factors 
is consistent with half-lives calculated in the spherical WKB or DWBA approximations. 
Recently reported proton decay rates indicate that the missing region of light rare-earth 
nuclei contains strongly deformed nuclei at the drip-lines. The lifetimes of deformed proton 
emitters provide direct information on the last occupied Nilsson configuration, and therefore 
on the shape of the nucleus. Modern models for proton decay rates from deformed nuclei 
have only recently been developed |§. However, even the most realistic calculations are not 
based on a fully microscopic and self-consistent description of proton unstable nuclei. In 
particular, such a description should also include important pairing correlations. 

The model that we use for the description of ground-state properties of proton emitters, 
is formulated in the framework of Relativistic Hartree Bogoliubov (RHB) theory. Models of 
quantum hadrodynamics that are based on the relativistic mean-field approximation, have 
been successfully applied in the description of a variety of nuclear structure phenomena in 
spherical and deformed /3-stable nuclei ||, and more recently in studies of exotic nuclei far 
from the valley of beta stability. RHB presents a relativistic extension of the Hartree- Fock- 
Bogoliubov theory. It was derived in Ref. in an attempt to develop a unified framework 
in which relativistic mean-field and pairing correlations could be described simultaneously. 
Such a unified and self-consistent formulation is especially important in applications to drip- 
line nuclei, in which the separation energy of the last nucleons can become extremely small, 
the Fermi level is found close to the particle continuum, and the lowest particle-hole or 
particle-particle modes are embedded in the continuum. The RHB model with finite range 
pairing interactions has been used to describe halo phenomena in light nuclei 0, the neutron 
drip-line in light nuclei |7]], the reduction of the spin-orbit potential in drip-line nuclei ||, 
ground-state properties of Ni and Sn isotopes ||, the deformations and shape coexistence 
phenomena that result from the suppression of the spherical N = 28 shell gap in neutron-rich 
nuclei [ID]. In particular, in Ref. ]TT[ we have applied the RHB model to describe properties 



of proton-rich spherical even-even nuclei 14 < Z < 28 and N = 18, 20, 22. It has been 



shown that the model correctly predicts the location of the proton drip-line. The isospin 
dependence of the effective spin-orbit potential has been studied. 

In the relativistic Hartree-Bogoliubov model, the ground state of a nucleus |$ > is rep- 
resented by the product of independent single-quasiparticle states. These states are eigen- 
vectors of the generalized single-nucleon Hamiltonian which contains two average potentials: 
the self-consistent mean-field T which encloses all the long range particle-hole (ph) corre- 
lations, and a pairing field A which sums up the particle-particle (pp) correlations. The 
single-quasiparticle equations result from the variation of the energy functional with respect 
to the hermitian density matrix p and the antisymmetric pairing tensor k. In the Hartree 
approximation for the self-consistent mean field, the relativistic Hartree-Bogoliubov (RHB) 
equations read 

,i„- m -x A \»)\ 

V "A* -h D + m + \J \V k (r)J \V k (r)J 
where hp is the single-nucleon Dirac Hamiltonian Q], and m is the nucleon mass 

h D = -ia ■ V + P(m + g a a{r)) + ^r 3 ^°(r) + g p p°(r) + e (1 ~ Ts) A°(r). (2) 

The chemical potential A has to be determined by the particle number subsidiary condition 
in order that the expectation value of the particle number operator in the ground state 
equals the number of nucleons. The column vectors denote the quasi-particle spinors and 
Ej. are the quasi-particle energies. The Dirac Hamiltonian contains the mean-field potentials 
of the isoscalar scalar cr-meson, the isoscalar vector a;-meson, the isovector vector p-meson, 
as well as the electrostatic potential. The RHB equations have to be solved self-consistently 
with potentials determined in the mean-field approximation from solutions of Klein-Gordon 
equations. The equation for the sigma meson contains the non-linear a self-interaction 
terms. Because of charge conservation only the third component of the isovector p-meson 
contributes. The source terms for the Klein-Gordon equations are calculated in the no-sea 
approximation. In the present version of the model we do not perform angular momentum 
or particle number projection. 
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The pairing field A in ([!]) is an integral operator with the kernel 

A a6 (r, r') = -^] V abcd (r, r')n cd (r, r'), (3) 

c,d 

where a, b, c, d denote quantum numbers that specify the Dirac indices of the spinors, 
V a b C( i(r,r') are matrix elements of a general relativistic two-body pairing interaction, and 
the pairing tensor is defined 

«w(r,r')= E U: k (r)V dk (r>). (4) 

E k >0 

In most of the applications of the RHB model we have used a phenomenological non rela- 
tivistic interaction in the pairing channel: the pairing part of the Gogny force 

V PP {1,2) = Yl e- ((ri - r2)/ ^ )2 (Wi + BiP a - HiP T - MiP a P T ), (5) 

i=l,2 



with the set D1S [|12j for the parameters /ij, Wi, B>i, Hi and M* (i = 1,2). This force has been 



very carefully adjusted to the pairing properties of finite nuclei all over the periodic table. 
In particular, the basic advantage of the Gogny force is the finite range, which automatically 
guarantees a proper cut-off in momentum space. 

The eigensolutions of Eq. ([[]) form a set of orthogonal and normalized single quasi- 
particle states. The corresponding eigenvalues are the single quasi-particle energies. The 
self-consistent iteration procedure is performed in the basis of quasi-particle states. A simple 
blocking procedure is used in the calculation of odd-proton and/or odd-neutron systems. The 
resulting quasi-particle eigenspectrum is then transformed into the canonical basis of single- 
particle states, in which the RHB ground-state takes the BCS form. The transformation 
determines the energies and occupation probabilities of the canonical states. 

The input parameters are the coupling constants and masses for the effective mean-field 
Lagrangian, and the effective interaction in the pairing channel. As in most applications of 
the RHB model, we use the NL3 effective interaction |13[ for the RMF Lagrangian. Proper- 
ties calculated with NL3 indicate that this is probably the best RMF effective interaction so 
far, both for nuclei at and away from the line of /^-stability. For the pairing field we employ 



the pairing part of the Gogny interaction, with the parameter set D1S [TJ 



In Fig. [| we display the one-proton separation energies 

S P (Z,N) = B(Z,N)-B(Z-1,N) (6) 

for the odd-Z nuclei 59 < Z < 69, as function of the number of neutrons. The model 
predicts the drip-line nuclei: 124 Pr, 129 Pm, 134 Eu, 139 Tb, 146 Ho, and 152 Tm. In heavy proton 
drip-line nuclei the potential energy barrier, which results from the superposition of the 
Coulomb and centrifugal potentials, is relatively high. For the proton decay to occur the 
odd valence proton must penetrate the potential barrier, and this process competes with 
(3 + decay. Since the half-life for proton decay is inversely proportional to the energy of the 
odd proton, in many nuclei the decay will not be observed immediately after the drip-line. 
Proton radioactivity is expected to dominate over f3 + decay only when the energy of the odd 
proton becomes relatively high. This is also a crucial point for the relativistic description of 
proton emitters, since the precise values of the separation energies depend on the isovector 
properties of the spin-orbit interaction. 

The calculated separation energies should be compared with recently reported experi- 
mental data on proton radioactivity from 131 Eu, 141 Ho ||, 145 Tm |TJ]], and 147 Tm [JOSfl . The 



131 Eu transition has an energy E p = 0.950(8) MeV and a half-life 26(6) ms, consistent with 
decay from either 3/2 + [411] or 5/2 + [413] Nilsson orbital. For 141 Ho the transition energy 
is E p = 1.169(8) MeV, and the half-life 4.2(4) ms is assigned to the decay of the 7/2" [523] 
orbital. The calculated proton separation energy, both for 131 Eu and 141 Ho, is of —0.9 MeV. 
In the RHB calculation for 131 Eu the odd proton occupies the 5/2 + [413] orbital, while the 
ground state of 141 Ho corresponds to the 7/2^ [523] proton orbital. This orbital is also oc- 
cupied by the odd proton in the calculated ground states of 145 Tm and 147 Tm. For the 
proton separation energies we obtain: —1.46 MeV in 145 Tm, and —0.96 MeV in 147 Tm. 
These are compared with the experimental values for transition energies: E p = 1.728(10) 
MeV in 145 Tm, and E p = 1.054(19) MeV in 147 Tm. When compared with spherical WKB 
or DWBA calculations [ 16] , the experimental half-lives for the two Tm isotopes are con- 
sistent with spectroscopic factors for decays from the hu/ 2 proton orbital. Though our 
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predicted ground-state configuration 7/2~[523] indeed originates from the spherical h u / 2 
orbital, we find that the two nuclei are deformed. 145 Tm has a prolate quadrupole deforma- 
tion fa = 0.23, and 147 Tm is oblate in the ground-state with (3% = —0.19. Calculations also 
predict possible proton emitters 136 Tb and 135 Tb with separation energies —0.90 MeV and 
— 1.15 MeV, respectively. In both isotopes the predicted ground-state proton configuration 
is 3/2+[411]. 

The calculated mass quadrupole deformation parameters for the odd-Z nuclei 59 < Z < 
69 at and beyond the drip line are shown in Fig. ^|. Pr, Pm, Eu and Tb isotopes are strongly 
prolate deformed (02 ~ 0.30 — 0.35). By increasing the number of neutrons, Ho and Tm 
display a transition from prolate to oblate shapes. The absolute values of j3 2 decrease as we 
approach the spherical solutions at N = 82. The quadrupole deformations calculated in the 
RHB model with the NL3 effective interaction, are found in excellent agreement with the 



predictions of the macroscopic-microscopic mass model |17 . 

A detailed analysis of single proton levels, including spectroscopic factors, can be per- 
formed in the canonical basis which results from the fully microscopic and self-consistent 
RHB calculations. For the Eu isotopes this is illustrated in Fig. |3], where we display the 
proton single-particle energies in the canonical basis as function of the neutron number. 
The thick dashed line denotes the position of the Fermi level. The proton energies are 
the diagonal matrix elements of the Dirac hamiltonian hp (||D in the canonical basis. The 
phase-space of positive-energy states should not be confused with the continuum of scat- 
tering states which asymptotically behave as plane waves. The RHB ground-state wave 
function can be written either in the quasiparticle basis as a product of independent quasi- 
particle states, or in the canonical basis as a highly correlated BCS-state. In the canonical 
basis nucleons occupy single-particle states. The canonical states are eigenstates of the RHB 
density matrix. The eigenvalues are the corresponding occupation numbers. In particular, 
we notice that for the proton emitter 131 Eu, the ground-state corresponds to the odd valence 
proton in the 5/2 + [413] orbital. 

In conclusion, this study presents a detailed analysis of deformed proton emitters 59 < 



Z < 69 in the framework of the relativistic Hartree-Bogoliubov theory. We have investigated 
the location of the proton drip-line, the separation energies for proton emitters beyond the 
drip-line, and ground-state quadrupole deformations. The NL3 effective interaction has 
been used for the mean-field Lagrangian, and pairing correlations have been described by 
the pairing part of the finite range Gogny interaction D1S. The RHB results for proton 
separation energies are found to be in very good agreement with recent experimental data 
on direct proton decay of 131 Eu, 141 Ho, and 147 Tm. The theoretical value is not so good 
for 145 Tm; the calculated and experimental proton energies differ by more than 200 keV. 
Predictions for the deformed single-particle orbitals occupied by the valence odd protons, 
are consistent with experimental half-lives for proton transitions. The model also predicts 
possible proton emitters 136 Tb and 135 T. The present analysis has shown that the RHB with 
finite range pairing provides a fully self-consistent microscopic model which can be used to 
map the entire proton drip-line for medium-heavy and heavy nuclei 51 < Z < 83. 
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Figure Captions 

FIG. 1. Calculated one-proton separation energies for odd-Z nuclei 59 < Z < 69 at and beyond 
the drip-line. 

FIG. 2. Self-consistent ground-state quadrupole deformations for the odd-Z nuclei 59 < Z < 69 
at the proton drip-line. 

FIG. 3. The proton single-particle levels for the Eu isotopes. The dashed line denotes the 
position of the Fermi level. The energies in the canonical basis correspond to ground-state solutions 
calculated with the NL3 effective force of the mean-field Lagrangian. The parameter set D1S is 
used for the finite range Gogny-type interaction in the pairing channel. 
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